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Edwards- Wilkinson type models are studied in 1 + 1 dimensions and the time-dependent distri- 
bution, Pl{vj'^ ,t), of the square of the width of an interface, , is calculated for systems of size L. 
We find that, using a flat interface as an initial condition, PL{w^,t) can be calculated exactly and it 
obeys scaling in the form {w'^)cx,Pl{w^ ,t) = <l>(w^/ {w'^)ao,t/L'^) where {w'^)ac is the stationary value 
of . For more complicated initial states, scaling is observed only in the large-time limit and the 
scaling function depends on the initial amplitude of the longest wavelength mode. The short-time 
limit is also interesting since Pl {w^ , t) is found to closely approximate the log-normal distribution. 
These results are confirmed by Monte Carlo simulations on a 'roof-top' model of surface evolution. 
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I. INTRODUCTION 

Interfaces play an important role in a number of phys- 
ical, chemical, and biological phenomena. Their fluctu- 
ations display universal features and, accordingly, mod- 
els of interface motion have been analyzed and distin- 
guished in terms of universality classes The clas- 
sification usually proceeds by taking systems of various 
sizes L and measuring the time-evolution of the average 
of the square width, (w^), of the interface. Then the val- 
ues of the static {Q and dynamic (z) exponents which 
determine the universality class are obtained ||^ by ob- 
serving collapse of data in accordance with the scaling 
form (^2) 

In practice, this procedure is not so easy to realize 
since {w'^) is an integral over all modes in the system and, 
consequently, large corrections to scaling are present. Re- 
cently, it was suggested Q that an alternative and more 
detailed characterization of interfaces may be obtained 
through the probability distribution of the random vari- 
able w"^. The steady-state distribution function pI^^\w'^) 
has been calculated exactly for various growth models 
[^j-^ and it has been found that P^'^'' {w^ ) defines a scal- 
ing function, 



(1) 



which is a universal characteristic of the interface fiuctu- 
ations. 

The above scaling function, however, distinguishes 
only among static universality classes. For example, 
$*^^-'(a;) is the same for both the one-dimensional {d — 1) 
Edwards- Wilkinson (EW) |^ and Kardar-Parisi-Zhang 
(KPZ) (t) equations since these models describe processes 
which have the same steady-state distributions and dif- 
fer only in the scaling of their dynamics. Thus, in order 
to distinguish among dynamical universality classes one 
should extend the static results to the time-dependent 



width distributions, PL{w'^,t). This is what we shall do 
here for the simplest model of surface growth, the d — 1 
EW equation 

Since the EW equation is linear, much of the calcula- 
tion can be done analytically (Section ||) and, in partic- 
ular, we can derive Pl{w^, t) and the associated dynam- 
ical scaling function (w-^) oo Pl (w^ ,t) in closed form for 
the case of a flat initial surface. An interesting feature 
of the result is that the short-time limit is closely ap- 
proximated by the log-normal distribution (Section III). 
Arbitrary initial conditions are harder to treat and we ob- 
tain general results only for the long-time limit where the 
dependence on the initial-state disappears except for the 
initial amplitude of the longest wavelength mode (Section 

In order to carry out a limited check of the univer- 
sality of our results we used Monte Carlo simulations to 
study a 'roof-top' model of surface evolution This 
model belongs to the EW universality class when the 
overall velocity of the surface is zero. Otherwise, it is in 
one universality class with the KPZ equation. Excellent 
agreement is found (Section [V ) between the analytic and 
simulation results for the dynamical scaling function in 
the EW limit and, furthermore, we also find that this 
scaling function is easily distinguishable from the corre- 
sponding function obtained for the KPZ case. 



II. CALCULATION OF THE WIDTH 
DISTRIBUTION 

A simple model of surface evolution governed by sur- 
face tension and noise is the EW equation M: 



dh{x, t) 
dt 



(2) 



Here h{x, t) is the height of the surface at sites < 2; < L, 
I' is a constant related to the dynamical surface tension. 
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and 7y is a Gaussian white noise of strength F: 

{ri{x, t)r]{x' ,t')) = 2rS{x - x')5{t -t') . (3) 

For simphcity, we shaU assume periodic boundary condi- 
tions. It should be noted, however, that free boundary 
conditions may be more reahstic in higher dimensions 
where comparisons with experiments are possible. 

Our aim is to calculate the time-dependent distribu- 
tion, Pi,{vlP' , t), for an arbitrary initial condition, h{x,Q). 
The derivation follows along the line that has been 
worked out for the static case ^ with extra complica- 
tions arising fr'om time dependence as well as from initial 
conditions. 

The quantity w^(t) is defined for a configuration h{x, t) 
as the mean square fluctuations of the height: 



(4) 



where the time-dependent average, f(t), of a fimction 
f{x,t) is obtained as its spatial average 



(5) 



/(t) = - / dxf{x,t) 



The first step of the calculation of Pl {'uP' , t) is writing it 
in terms of a path integral 

PL{w'',t)= I V[h] Slw"" -(p-t)]p{{h},t) (6) 



where p({/i},t) is the path-probability that the surface 
evolves from an initial state h{x, 0) to a configuration 
t) in time t. The dependence on the initial condition 
is not written explicitly though it is understood that this 
dependence is important for any finite t. 

The Laplace transform of equation ^ gives the gen- 
erating function of the moments of P^ (w^ , t) : 

GL{X,t)^ / dCPL{C,t)e-^^ , (7) 
Jo 

and one finds that G(A, t) is the following path integral 
GL(A,t)= [ V[h] p{{h},t)exp\-X{h^ -h'^)] . (8) 



The next step is to note that the above path integral 
can be written as an infinite product of ordinary integrals 
provided the system is described in terms of Fourier am- 
plitudes. Indeed, let us write 



h{x,t) - h{t) 



E 



Cn{t)e 



iknX 



(9) 



where fc„ = 27rn/L and c_„ = c* (note that cq = 0, 
thus the n — Fourier mode can be left out from further 
considerations). The EW equation is replaced now by an 
infinite set of ordinary differential equations: 



(10) 



where 77,1 (i), defined as the Fourier transform of ■i]{x,t), 
is also an uncorrelated white noise. We can see that the 
Fourier modes with different n-s are decoupled and evolve 
independently. Thus the probability of a path that an 
initial state characterized by a set of {c„(0)} evolve into 
{cnit)} is just the product of probabilities Pn[cn(t)\cn{0)] 
that c„(0) evolves into c„(t): 



p[{c„(t)}|{c„(0)}] = n PnKit)\CniO)] 



(11) 



n— — 00 



Equation (|T^) is well known as the Langevin equation 
of Brownian motion and thus Pn[c„(t)|c„(0)] can be ob- 
tained from the Fokker-Planck description of the process 



{t)\cnm 



where: 



exp 



(c„(<)) =c„(0)e-'"="' , niit) 



r(i 



-2vk'it\ 



Lvk3, 



(12) 



(13) 



Now we can write the functional integral (||) as a prod- 
uct of integrals over the coefficients c„: 



GL{\t) = n y dcJclpl [c„(t)|c„(0)]e-2^l' 



(14) 



Substituting p[c„(t)|c„(0)] from (p^), the Gaussian inte- 
grals can be calculated and the inverse Laplace transform 
of the result gives the time-dependent probability distri- 
bution in a scaled form which does not depend explicitely 
on L: 



{w^)ooPL{v?,t) EE $(x,T,{s„o}) = 



dy 
27rj 



n 

n=l 



exp 



-y 4o e-" 7(1 + 2/a„) 



{1 + yan) 



(15) 



Here {w'^)^= LT /{I2v) and 
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r(l 



(ttti)^ 

with the scaling variables given by 



jJ2 ' nO 



,2 _ 



2|c„(0)|^ 



(16) 



(17) 



In case of flat initial surface (s„o = 0), one can evaluate 
the integral ( p^ exactly by collecting contributions from 
simple poles at — l/a„ and one finds a scaling function of 
two variables: 



2 



^{x,t) = $(X,T,{0}) = 

? n 



(18) 



For any finite x, the above sum can be approximated by 
a finite number of terms and $(a:,r) can be evaluated 
with a given accuracy. Difficulties arise only in the limit 
of r ^ where $(x,t) goes into a delta-function. The 
evolution of $(a;, r) obtained by numerical evaluation of 
( p^ ) can be seen in Fig. |^. 
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FIG. 1. Evolution of the scaling function of width distri- 
bution in case of a flat initial interface [eq.(I8)]. The scaling 
variables x and r are given in eq.(17). Analytical results for 
the EW model are compared with log-normal fits and with 
Monte Carlo results on the 'roof-top' model. The exact results 
and the lognormal fits are indistinguishable within linewidth 
for T < 4ri. 

We conclude this section by calculating the time de- 
pendence of (w^) which will be needed for comparing the 
time-scale of MC simulations with the time-scale of the 
EW equation. We find 
tion as — 9aG'(A)|a=o : 



(sIq Sioo) initial values for s^o ^^d then {w'^)t ap- 
proaches its steady-state value from below and above, 
correspondingly. If the two values converged, one may 
assume that the steady state has been reached. This 
type of checking for steady state is widely used in sim- 
ulations of equilibrium systems such as the Ising model 
where completely ordered and disordered initial states 
are employed. Similar procedures, however, do not seem 
to have been followed in the simulations of surface evo- 
lution models. 

For short times, one can change the sum (|9|) into an 
integral and {'w'^)t — {w^)t=o 



=n can also be calculated 



V 



(21) 



The above result is valid for a flat initial state as well as 
for an initial surface containing finite number of nonzero 
Fourier terms (s„o 7^ for n < Umax)- 



III. LONG- AND SHORT-TIME ASYMPTOTICS 

The case of arbitrary initial conditions with nonzero 
Sno-s is complicated by the presence of essential singu- 
larities in the function which is integrated in (p^. As a 
consequence, we are not able to evaluate the probability 
distribution in general. The representation ( p^ ) is useful, 
however for finding the long- and short-time asymptotics 
of <i>(a;, t). 

For r ^ 00, one has to keep contributions which are 
proportional to e""^ ^ 1. It follows then that $ depends 
on the initial state only through the initial amplitude, 
sio, of the longest wavelength mode. The calculation of 
the terms which are proportional to e""^ involves collect- 
ing contributions from both simple poles and quadratic 
singularities in ( p^ ) with the result that the siq depen- 
dence appears in a prefactor in front of a new scaling 
function 



[w'^)t from the generating func- ^{x,T,{sno}) — 



= 1 



00 

E 

ri=l 



o2 ^2 



1 - 



^(x) e-^ + 0{e-'^) , (22) 



(19) 



where 



In the long time limit (e <C 1) this expression reduces 
to 



mix) — — X exp < x 

^ ' 3 U 6 / ^1 6 



l + {4^-±)e-^ + 0{e-^n , (20) 



9 00 A 
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m=2 



exp — —m x 



(23) 



and one can see that {w'^)t approaches its steady-state 
value much faster if we start with an initial surface where 



sJq is set to its steady-state value sfg 



*loo 



6/7r2 



The above equation also suggest a method for finding 
out in simulations if the system has settled to its steady- 
state. One can choose small (s^g ^ sf^) and large 



This scaling function ^'(a:) is shown in Fig. ^ and dis- 
cussed in Section [V, 

Comparing equation (^) with (pO|), one can see that 
the same prefactor, s^q — G/tt^, appears in front of e~'^ in 
both cases. Thus the relaxation of the scaling function 
is also accelerated if siq is chosen to be the steady-state 
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value sioo- Furthermore, it also follows that the steady- 
state distribution can also be bracketed by choosing small 
and large initial values for sio- 

Now we turn to the description of the short-time limit 
of <i>(a;,T) in case of a flat initial condition. The de- 



scription is based on an earlier observation |11| that the 
fluctuations of chemical reaction fronts which are sup- 
posed to belong to the EW universality class produce a 
$(x,t) which is rather well approximated by log-normal 
distribution 



$(a;,T) « C{x,xo,(j) = 



1 r \n'{x/xo) \ 




(24) 



where xo{t) and a{T) are fitting parameters which can be 
determined from various considerations. We have deter- 
mined xo{t) and (t(t) by equating both the maxima and 
positions of the maxima Xm of the two functions t) 
and C(x, XQ,a). The values of x^ and t) — <i>„j(T) 

have been determined numerically from (|l8|) and then a 
was obtained by solving the following equation: 



V2^x.m^„,a = e 



(25) 



and finally, Xq was expressed as Xo = Xm,exp{a'^). For 
sufficiently short times, the width of the distribution goes 
to zero and cr — > 0. In this limit the expressions for xq 
and a simplify to xq « 



and a w l/{y/2TrXm^m) 
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FIG. 2. The analytical results and the log-normal fits of 
Fig.l on a log-log plot. 

The results of this fitting procedure can be seen on 
Fig. |l|. Qualitatively, the fit is quite good over the whole 
time interval and it clearly becomes excellent at short 
times. In order to make the quality of the short-time 
fit more apparent, we have redrawn the curves of Fig. |l| 
on a log-log plot with the results shown on Fig. ^. As 
one can see again, the log-normal fit becomes better at 



small times and, at r = 0.039, the fit becomes practi- 
cally indistinguishable from ^(x,t) in an interval where 
the function decreases from its maximum value by three 
orders of magnitude. 

A more quantitative description of the quality of the 
log-normal fit can be given by defining a relative distance 
between the two functions as 



e{r) 



|$(g;,r) - C{x,XQ,a)\ 

$m(T) 



(26) 



This distance increases with r and reaches its maximum 
value Imax = ^(oo) « 0.14 in the stationary state. For 
T -> 0, we find that £ - r^/^ and £ < 0.01 for r < 
0.02. Although the diminishing relative distance actually 
comes from the ratio of a strongly divergent <I>„j(t) ^ 
T~^/^ and a less divergent maximum distance maxj; |(E> — 
C\ ^ r^^/"*, plots of the two functions which extend from 
zero to their maxima are indistinguishable (see Fig.Q) for 
T < 0.1. 



IV. MONTE CARLO SIMULATIONS 

In order to see if the dynamic scaling found for the 
width distribution had the expected universality, we car- 
ried out Monte Carlo simulations for a 'roof-top' model 
of surface evolution [^,^. In this model the height of 
the surface is characterized by a single-valued function 
hi at sites i = 1,2,...,L and periodic boundary condi- 
tions hi-i-L = hi are imposed. The height differences are 
restricted to /i^+i — hi = ±1 and the evolution consists 
of particles being deposited at local minima or evaporat- 
ing from local maxima of the surface with rates p-^- and 
P- = 1 — p+, respectively. If J3+ = P- — 1/2, the model 
belongs to the universality class of EW model while for 
p+ 7^ p_ the universality class is that of the KPZ equa- 
tion 1^. 

For equal rates, one can obtain |^ an exact expression 
for {'w'^)t and comparing the result with the solution of 
the EW equation (|9|), the time-scale of the MC simu- 
lation can be related to that of the EW equation. In 
this way one finds that the parameters and F should 
be set to = r = 1/2. Then the x-s and r-s in the 
MC data and in the EW equation are related in a unique 
way and there is no parameters to fit when the $-s are 
compared. Fig. |l| shows both the ^(x, t)-s obtained from 
simulation and the theoretical curves of the EW model. 
One finds good agreement although a small systematic 
shift of the MC curves towards larger values of x can 
be observed. This shift is due to the fact that, in the 
'roof-top' model, the initial surface is not entirely flat 
(w^t^o = 1/4) in variance with the w^^q = used in the 
theoretical calculation. This difference should disappear 
in the L — > oo limit and, indeed, one can see that the 
difference is smaller for the L = 512 sample as compared 
to the L = 256 system. 
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We have also examined the function ^{x) which char- 
acterizes the scahng of the long time relaxation of the dis- 
tribution function in the EW model (|2^). Since one can 
find a large enough time- window where {'uP')r — {w'^)oc ^ 
e""^ for all x, the function "^{x) can be determined ac- 
curately and, as can be seen from Fig. there is an 
excellent agreement with the theoretical curve. Thus we 
can conclude that the 'roof-top' model which belongs to 
the EW universality class for p+ = = 1/2 indeed pro- 
duces the same time dependent distribution $(a;,r) as 
the EW equation. 
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this moment we have results only for the one-dimensional 
EW and KPZ processes but there does not seem to be 
any principal difficulty in extending these calculations to 
other processes and to higher dimensions by using exact 
solutions, renormalization-group methods, and simula- 
tions. 

An interesting byproduct of our calculation is the re- 
sult that the early-time width distribution in the EW pro- 
cess is practically identical to the lognormal distribution. 
Lognormal-like distributions tend to emerge more often 
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in biological and social sciences than in physics [p2 13 
and they are usually understood in terms of the 'law of 
proportionate effect' or on the basis of the assumption 
that an event occurs only if a large number of indepen- 
dent 'sub-events' take place. In our case, the lognormal 
distribution is produced by EW dynamics and it appears 
as a characteristic of the initial roughening of an inter- 
face. Whether this generation of lognormal-like distribu- 
tions was new or it was equivalent to one of the standard 
derivations remains to be understood. 
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FIG. 3. Scaling function, '^'{x), describing the long-time 
relaxation of the width distribution [eq.(22)]. The theoretical 
curve for the EW model is compared with Monte Carlo results 
for both the EW and KPZ limits of the 'roof-top' model. The 
initial states are either flat or contain a single sine perturba- 
tion with the longest available wavelength. 

In order to investigate if the ^(x, t)-s characterizing 
the EW and the KPZ classes were distinguishable, we 
have also studied the long time behavior of the 'roof- 
top' model for unequal rates (p+ = 1, p_ = 0). In this 
case, we find that (w^) oo PL{x,t) — $(x,oo) decays with 
time exponentially, exp(— a^t), with a relaxational rate, 
aL independent of x and, as can be seen on Fig. ^ the 
coefficient, ^'(a;), of the exponential differs significantly 
from that of the corresponding EW scaling function. 



V. FINAL REMARKS 

It has been demonstrated previously |3|-||] that one can 
build a 'picture gallery' of scaling functions for steady- 
state width distributions and this gallery may be used 
for distinguishing the static universality classes of growth 
processes. Here we have made the first steps towards 
building a similar gallery for dynamic scaling functions 
and we believe that this gallery will be equally instru- 
mental in recognizing dynamical universality classes. At 
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